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Abstract 



We propose a simple idea to construct ID integrable models with impu- 
rities. We illustrate the strategy for a supersymmetric t — J Hamiltonian 
in considerable detail. The impurity comprises the local deformation of the 
hopping and exchange integrals as well as a three-body charge-current inter- 
actions on neighboring sites. We explore the thermodynamic properties of 
the system at low and high temperatures, and obtain results that are beyond 

boundary conformal-field theory. 
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Recently there has been new developments in the theory of impurities in Fermi and Lut- 
tinger liquids. This renewal of interest in the quantum impurity problems has been stim- 
ulated in part by the search for non-Fermi-liquid fixed-points beyond one space dimension 
and its plausible relevance to the theory of heavy fermions and high temperature supercon- 
ductors. In this context, Affleck, Ludwig and collaborators have used the (multi-channel) 
Kondo problem as a laboratory to test, develop and extend significantly the techniques of 
conformal field theory in the presence of boundaries [|l],^. These powerful methods allow 
to obtain single-particle as well as thermodynamic properties near the critical point. Using 
simpler bosonization and renormalization group techniques, Kane and Fisher have shown 
that a potential scatterer embeded in a Luttinger liquid is driven to a strong-coupling fixed 
point by the electron-electron repulsive interactions 0. Emery and Kivelson were able 
to map the two channel Kondo problem to a resonant-level Hamiltonian which could be 
solved exactly at a particular point Fabrizio and Gogolin have generalized the Emery- 
Kivelson solution to the isotropic four channel and together with Nozieres to the anisotropic 
multi-channel Kondo problem [§,0. Combining bosonization and conformal field theory 
results, Sengupta and Georges have proposed a solution for a two-channel two-impurity 
Kondo Hamiltonian [0. Significant progress in a different direction has also been achieved 
by Ghoshal and Zamolodchikov. These authors were able to solve exactly integrable field 
theories with boundaries Interesting applications of these results have been performed 
by Fendley et al. and Tsvelick 0. 

Despite this important progress, the problem of an (few) impurity(ies) (potential, mag- 
netic or boths) embeded in a strongly correlated system is not well understood. There is 
a notable exception: a spin S > 1/2 impurity in a spin 1/2 Heisenberg chain solved many 



years ago by Andrei and Johannesson |T0| and generalized to arbitrary spins by Lee and 



Schlottmann, and Schlottmann [|TT| . Very recently, it was argued in reference |]T2[ that these 
models correspond to non-generic fixed points. In the present letter, we want to explore a 
new idea that allows us to construct exactly solvable models with a (few) defect (s). Short 
range interactions between the impurities can also be incorporated. There is a prize to pay, 



however, namely some terms in the Hamiltonian are difficult to justify on physical grounds. 
Their presence is required by formal reasons (integrability) and should be generically irrel- 
evant to the low-energy properties of the system. The idea is simple and can be explained 
as follows. Consider an integrable model that can be formulated within the framework of 
the Quantum Inverse Scattering Method (QISM) |]T1|. In the QISM, the Hamiltonian of the 
model is written as the logarithmic derivative of an homogeneous transfer matrix at a special 
value of the spectral parameter (an infinite number of conserved kinematic quantities can be 
generated by performing successive derivatives of the logarithm of the transfer matrix). In 
general, however, the transfer matrix may depend on a set of (unequal) parameters. In the 
statictical-mechanics language this corresponds to an inhomogeneous vertex model. Define 
the Hamiltonian as the logarithmic derivative of the corresponding inhomogeneous transfer 
matrix. Provided we restrict to an (few) inhomogeneity(ies), we have a model with a (few) 
defect (s) that is exactly integrable by construction. We have performed preliminary calcula- 
tions in a variety of models: spinless fermions, spin chains and strongly correlated systems. 
A general discussion will appear elsewhere ||T^|. 



In this work, we want to illustrate the idea by focusing on one particular model that 
is a prototype for strongly correlated electrons, i.e., the t — J model. In the absence of 
impurity, the t — J model is known to be solvable at the supersymmetric point (away from 
the integrable point, much is also known through numerical studies, see ref. [ff^) and various 



properties of the model have been investigated in the recent past |T7|-[19|. To construct at — J 
model with impurity, we follow the strategy proposed above and arrive at the Hamiltonian: 
H = Hq + Hj, where Hq = Yln=i Hnn-i denotes the usual t — J Hamiltonian, 

Hnn+l = Vi- Y,{ciaCn+la + h.C.) + J(S„+i ■ S„ - ^"^"'^^ ) + + nn+l}V , (1) 

where V projects out the doubly occupied sites, fi„ = X^o- cl^a^na, J = 2 and periodic 
boundary conditions have been assumed. The impurity contribution can be written as 



— Im+l ■^mm—l -^mm+l — lm+1 ; -^mm— l] 7 (2) 
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where the couphng constant < g < 1 (this restriction follows from the hermiticity of H) 
and m denotes the center of the defect that involves, besides m, its two adjacent sites m± 1. 
The first term on the r.h.s. of (^) enhances antiferromagnetic correlations between mil 
sites while ferromagnetic alignment is favored between m and its two adjacent sites. The 
second term corresponds to a three-body interaction, 

^ {(1 - ^m) EOr+lm-1 + h.C.) - (1 - EOr+l™ + h.C.) - (1 - Um+l) + h.C.) 

a cr a 

+ '^Z^^ma^ma' j m+lm-1 + h.C.) — '"^{cl^-lcrCm-la' jm+lm + h.C.) — ''^ic]^_^_i^Cm+la' jmm-1 + h.C.) 
era' era' aa' 

~'^i'^rna'^ma'Jm~-lm+l ~^ h.c)}V , (3) 
aa' 

with jX' = icla^ba' and J^' = ^ EA(cL'CaA)(4AC6a)- In (D, we have three types of terms: 

(i) a local charge current flowing between the sites m, m ± 1 provided one of them is empty; 

(ii) assisted spin-flip currents; (iii) pure assisted spin-currents. At g = 1, the three body 
term is absent and the model behaves as a chain of L — 1 sites with no impurity present. 

The eigenvalue problem associated to H is solved by the graded version of the QISM 
(see ref. [^). For convenience, we choose the Sutherland representation (see discussion in 
P^). The formal solution is casted into the form of the so-called Bethe ansatz equations. 



LpaiX]) + 5a(A,") = 2vr/; -EE <PU^1 - , (4) 

(3=c,s k=l 

where Pai^'j), for a = c,s, denotes the bare momentum with Ps{x) = 2tan^^(2x) and 
Pc{x) = 0, (paisi^) is the bare phase shift for the scattering of the a with the /? pseudoparticles 
(0°^(x) = — 2tan^^(x), 0°c(x) = (pcs{x) = 2 tan^^(2x),0°^(x) = 0) while 6a{x) denotes the 
bare phase shift for the scattering of the a pseudoparticle with the impurity, i.e., Ss{x) = 
2[tan-i(2[x-6'])-tan-i(2s)] and 4(x) = 0. The parameter 9 is defined by 6* = (g/il-g))^/"^. 
As usual the quantum numbers /" are integers or half-odd integers depending on the number 
of sites L, the number of particles N and the magnetization M in the system. The ground- 
state (of finite momentum ) and part of the low-energy spectrum are parametrized by real 



roots in equation (m (see ref. |1 
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Since we are interested in thermodynamic properties, we shall directly work at finite 
temperature with the thermodynamic potential per unit length H, T), where /i denotes 
the chemical potential, H the external magnetic field and T the temperature. At finite T, 
we have to consider all possible solutions of the set of equations (^). It turns out that the 
rapidities associated to the s-pseudoparticles cluster into strings in much the same way as 
they do in the Heisenberg chain while the "charge "rapidities remain real. Following 



Yang and Yang [^, we write down the thermodynamic potential and minimize with respect 
to the distribution of roots of the set of equations After some algebra (see for example 
I), the thermodynamic potential per unit length, i.e., the pressure, is found, up to an 



irrelevant constant, 

log[l-/(e„,,(x))] , (5) 



I f 



oo ^ r ]^ 

an{x) + — (a„(x -9) - a„(x)) 



where h = fisH, fiB is the Bohr magneton, P = (ksT)'^, a„(a;) = {27T)^^dpn{x)/dx with 
Pn{x) = 2 tan^^(2a;/'n,). Q naturally separates into a bulk contribution of 0(1) and an 
impurity contribution of 0(1/ L). The function f{esn) = (1 + e'^'^"")"^ coincides with the 
Fermi distribution function in the low-temperature limit while it differs from the former at 
finite T where the quasi-energies ens{x) are T-dependent. The latter are determined from 
the integral equations 

Pesn{x) = s*log(/(e 

sn—l )/( )) i^) 

+ SnA[s^\ogil- fie,))ix)-27iPJsix)] , (6a) 
Pecix)=Pfi-ph + J2an^\og{l-f{esn))ix) , (6b) 

n 

where s{x) = 1/2 cosh(7rx), -k denotes a convolution and the boundary conditions, 
limn^oo{^ns/n) = h, have been imposed. Notice that (^) and (|6ED do not depend on the 
impurity because the contribution of the latter to the free energy vanishes in the variation 
with respect to the densities of roots. 

In general, the set of coupled equations ( pa| ) and ( pbD are difficult to solve analytically. 



In the following, we solve these equations in the low- and high-T limits. We use techniques 



similar to those developped by Takahashi [|2^. We first discuss the high temperature region 
{(3J ^ 1), i.e., the weak couphng hmit. To lowest order in /5J, the pseudo-Fermi distribution 
functions /(cc) and /(e^n) are constants, and ( pa| ) and (|6BD simplify to /(csn) = F'"^ and 
/(ec) = 1 - where F„ = - - h'^) with z = e"^'* and 6 = [(1 + 



)/{l + z ^e^^)]^/^. Substitution of this result into the thermodynamic potential (|^) leads 



to: VIq = —f3^^ log [1 + 2e^^ cosh(/5/i) j . We note that the impurity does not contribute in 
this limit. Qq is just the thermodynamic potential per unit length for a lattice where each 
site has three possible states fl^. In this case the entropy, specific heat, etc. are easy to 
interpret. To higher order in j3, the quasi-energies are no longer constants and we expand 
them as Pe.jix) = Pe%ix) +j:n=ii^PJr6f (x), and Pe^ix) = Pe^ix) + j:n=ii^PJn^^Kx), 
where /3e°„ = log(F„_iF„_|_i) and /5e° = log^F^ / {1 — F^)) . Below we discuss the lightly doped 
Mott-Hubbard insulator. Near half-filling, we have solved for CiJH^) and Cj}\x) and their 



Fourier transforms are Cin (^) = ^{Fq/Fi] 



and 



(F„/F„„i)e-"l^ - (F„/F„+i)e-("+2)^ 
Q^\tu) = — )e^l'^l . Substituing this into (|), leads to the thermodynamic potential 
n = no + ni + 0{{f3J)^), where ^1 = ^1^ + n\/L = -{J/F^) (1 - g/L). From the latter we 
infer immediately that 

AS^ AC^} AM^ Ax^s ^Xc ^ ' 
where AS^''^ , AC*''^ , AM*''' , Ax^''^ ; ^Xc'* denote the 0{PJ) corrections to the impurity 
(host) entropy, specific heat, magnetization, spin susceptibility, and compressibility, respec- 
tively. The simplicity of this formula is remarkable. It expresses the fact that in the weak 
coupling limit, the impurity "mimics" the host chain. Notice that at g = 1, the expression 
for the thermodynamic potential density ( remember to multiply by L/{L — 1)) reduces to 
that of a chain of L — 1 sites. We have evaluated explicitly the entropy, specific heat, etc. to 
0{(3J) as a function of /i, h and T. For brevity, we mention below only the simplest results 
|T5| . Inverting the equation that defines the density p = 1 — 5, (5 = hole density) as function 
of p, yields for 5 ^ 1 



log 5 + log(l + 5) + log(2 cosh(/5/i)) + /3 J ( 1 - 1 



+ 0{5\{(3Jf 
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In the Sutherland representation, the Heisenberg chain is obtained for — oo ( remember 
that fi plays the role of the chemical potential for the holes). Substituing (j^) back into the 
entropy per unit length, at h = 0, we have 

-^^log2 + 5(l-[log5 + log2])-^ 
kB 2 

+ ^(1 - f ){5[log5 + log 2] + - 2[log5 + log 2])} + 0(5^ (pjf) , (9) 

The various contributions to the entropy are easy to interpret: (i) the entropy of free spins 
1/2 (of density 1 — 6); (ii) the entropy of free holes; (iii) a reduction of the entropy due to 
the exchange correlations (J > and log 5 < — log 2 for 5^1); (iv) an increase in entropy 
due to the impurity that depletes locally the exchange. The specific heat in zero field near 
half-filling behaves as 



C, = kBS[log6 + log2] 



[log 5 + log 2] - /3 J ( 1 - I) + ^ [log 5 + log 2] 



(10) 



The magnetization is finite in an external field and near half-filling takes the form 

:i - 6) tanh(/5/i) - pj (^1 - (1 - tanh2(/5/i))(l - 35 + 36^) + 0{6\ {(3Jf) , (11) 



M 



i.e., the magnetization is suppressed by (i) the doping that reduces the number of spins; (ii) 
the antiferromagnetic correlations due to the exchange. To next order, the defect enhances 
the magnetization locally, while the doping renormalizes down the effective exchange between 
nearest neighbors. The 5^ term can be interpreted as an attractive interaction between the 
holes that effectively increases the renormalized exchange and thereby penalizes the magnetic 
state. We have also evaluated the charge and spin susceptibilities to 0(5^). The general 
expressions are rather complex and will be discussed elsewhere [ p!5[| . 

We now pass to the low temperature region, i.e., strong coupling regime. We note that 



the quasi-energies es„(x) > for n > 2, h > (compare with ref. |2^) and so as T — > 



only esi(x) = es{x) and ec(x) survive in equations ( |5a] ) and ( |5b| ) [^. At T = 0, the integral 



equations determining es{x) and e^x) coincide with those of ref. |]18|. We have explicitly 



evaluated the distribution of roots as well as the quasi-energies in zero field near half-filling 



151 . The first non-trivial corrections to the spinon and holon quasi-energies appear to third 



order in the density of holes 60 {6 = 6q + 6i/L), i.e., 



ttJ 



2 cosh(7ra;) 



^0^ 



TT 



'C(3) 57rlog2 



2 log 2 



TT 
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[1 - 2tanh^(7ra;) 



(12a) 



TT 



Ps{x) ~ ^ - tanh"^(sinh(7ra;))(l + 60) 



+ 1 



TT 



3 V21og2 



TT log 2 



tanh(7ra;) 3C(3) / 2 log 2^ 
cosh(7ra;) 27r \ ■''^ 



(12b) 



eAx 



ttJ 



log 2 



TT 



- Rix) + ^ + 5? 



3C(3) , ,3 (3C(3) 




TT 



C(3)^,^, 5 log 2 ^„ 



2ti 12 log 2 



-R{x) - 



167r 



R'ix) 



PJx) ^ 27rF(x) 



TT 



27T 12 log 2 



^0^ 



TT 



TT 



3 12 log 2 



R'ix) , 



where Pc(s){x) denotes the momentum of the c(s)-type excitation relative to the ground- 
state, C(3) is the Riemann zeta function, R{x) = JJe [^'(l + ix/2) - *([1 + ix]/2)] /2'rr with 
the definition \l'(a:;) = d\ogr{x)/dx. R'{x) and R"{x) denotes the first and second deriva- 
tives of R{x) while 27iF{x) = (log7r)/2 + log (1(1 + ix/2)/r([l + ix]/2)). The 0(5'^) term 
in ec(x) is due to a chemical potential correction, i.e., /x = J[log2 — 3({3)Sq/4] + 0(5"^), 
that shifts the bottom of the holon spectrum upwards. From (|12a|) and ( |13a|) , we see that 
the spinous and holons behave as well defined objects with respective renormalized ve- 
locities Vc ~ 37r((3) J5o/8 log^ 2 and Vs ~ vrJ(l — 60) /2, to 0(5^). Thus, the holons and 
spinous scatter off the impurity one by one. This will be important when we deal with low- 
temperature transport properties []T3[. In zero field, the hole-charge is depleted in the region 
surrounding the defect, i.e., the impurity is screened by an amount 6i = —G{9)6o, where 
G{6) = 1 — T[R{9)/log2 > 0, while the magnetization vanishes locally, i.e., the impurity is 
magnetically inert. This redistribution of charge around the defect explains why the spinon- 
and holon-impurity phase shifts can simply be expressed in terms of the holon-holon and 



(13a) 
(13b) 



spinon-holon phase shifts [15 . 
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The zero-field ratios of impurity to liost susceptibilities can be shown to be xL/Xa = 
[P^ai.'ia) + Po(~5'o)]/2Pa('?a), i-G., cqual to the ratio of the average (over the left and 
right Fermi-points) density of states of the impurity to the host chain. Near half- 
filling, we find x\/x^ ~ ~^oG{6), i.e., the local magnetic response in zero field is con- 
trolled by the screening process, while for the compressibility ratio we have Xc/Xc ~ 
G{e) + (52(7r/21og2)2[(7r/21og2)i?"(^) + (37rC(3)/2 log^ 2)i?(^) + 2R\e)]. At finite magnetic 
field, the defect becomes magnetically actif. For weak fields, h <^ J, and near half- 
filling, the Wiener-Hopf method leads to a local magnetization that is linear in h, and 
Mi/Mh ^ e""^ -1- 5o[e^^ - (7ri?(^)/ log2)], for 9 < | log([e/27r]/i)|/7r. Notice that the 
impurity-bulk coupling is infrared finite in contrast to the Kondo problem (compare with 



ref. [^). At zero temperature, we can derive simple relations between the h = suscepti- 
bility ratios Xa/Xa their finite field counterparts, XaW/XaW, i.e.. 



X\{h) 



Xc Xs zi 



Vc f2Zs 
1 + — — - 
v., V2z., 



"1 r 



x; x^Vc 



Xc 



2Z. 



XsVs \2Zsc 



(14a) 
(14b) 



where Va, a = c,s, denotes the velocity of the a pseudoparticle at the T = Fermi points 
for finite h and the "dressed charge matrix" (essentially the S-matrix of the s- and 
c- pseudoparticles at the Fermi surface, see ref. [0). The formulas ( |14a| ) and ( |14b| ) are 
universal in the sense that only the coefficients in front of the matrix elements Za/B depend 
on the model under consideration ITS . 



We have evaluated the low-temperature corrections to the thermodynamic potential to 
0{T^), from which we have inferred the low-T magnetization, entropy, susceptibilities, a.s.o. 
We quote the specific heat = C^^^ + Cf^) + 0{T^), with 



CP 



' ' iPai(g°) + Pa.(-g°)~ 



1 Pm(g°) + Pm(-g°) 



(15a) 



^ L 2p„,(g0) 

pM) \ L 2p',M) y 



The ratios of zeroth, first and second derivatives of the Fermi surface density of states of 
the impurity to the host chain occur in (15a) and (15b). The coefficients Aa,Ba and Ca 



depend in a universal way on the Fermi surface velocity, the curvature of the spectrum and 
its rate of change as well as on the Scattering matrix and its higher derivatives. For the 
sake of brevity, we omit here these lengthy and complicated expressions ||15|. The linear 



term in T, (|15a]) , coincides with that inferred from the zero temperature finite size spectrum 
calculations [|l^], and can alternatively be derived from boundary conformal field theory P|. 
The next order term ( |15b|) cannot be calculated within the framework of conformal field 
theory. Very close to half-filling, equations ( |15a| ) and ( |15b| ) do not hold, however. At low 



T, when the charge carrier density 6o is small, a crossover takes place to a regime where the 
holon-liquid behaves as a non-degenerate ID Fermi system coupled to an impurity. A crude 
estimate for the onset of the crossover is obtained when the width of the holon "Fermi-sea" 
is of order fc^T, i.e., 5^ ^ [/c^T/J] (47r-i log 2 (SCiS))-^/'^). For instance, at T = lOK with 
an exchange of order 1 eV, 6c ~ 10^^ while at room temperature we have Sc ~ 10^^ |T5[. A 



similar crossover takes place for the spinon liquid near the transition to the ferromagnetic 
state. Following ref. I|,|l3 we can calculate the conductance through the impurity at low-T. 
We shall report on this elsewhere |15 . 



In summary, we have proposed a novel strategy to construct solvable models with impu- 
rities. The hamiltonians generated by our method contain terms that are difficult to justify 
on physical grounds. Yet some of those terms are, in the models considered so far, and 
should generically be, irrelevant to the low-energy properties of the model. We have illus- 
trated the idea by performing explicit calculations in a simple model with a single impurity. 
We have derived various results that are beyond boundary conformal field theory by solving 
the thermodynamic Bethe ansatz equations in the high-and low-temperature limits. The 
transport properties of these class of models will be explored in future work. 
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